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IS THE BURSTING RADIO-SOURCE GCRT J1745-3009 A DOUBLE NEUTRON STAR BINARY ? 
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ABSTRACT 

GCRT J1745-3009 is a peculiar transient radio-source in the direction of the Galactic Center. It 
was observed to emit a series of ~ 1 Jy bursts at 0.33 GHz, with typical duration ~ 10 min and at 
apparently regular intervals of ~ 77 min. If the source is indeed at the distance of the Galactic Center 
as it seems likely, we show that its observational properties are compatible with those expected from a 
double neutron star binary, similar to the double pulsar system J0737-3039. In the picture we propose 
the (coherent) radio emission comes from the shock originating in the interaction of the wind of the more 
energetic pulsar with the magnetosphere of the companion. The observed modulation of the radio signal 
is the consequence of an eccentric orbit, along which the separation between the two stars varies. This 
cyclically drives the shock inside the light cylinder radius of the less energetic pulsar. 

Subject headings: Radio continuum: stars — Radiation mechanisms: non-thermal — Stars: neutron 
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1. INTRODUCTION 

Recently Hyman et al. (2005) announced the remark- 
able discovery of a powerful bursting radio-source located 
about ldeg south of the Galactic Center (GC). GCRT 
J 1745-3009 was detected in 2002 during a radio transient 
monitoring of the GC region at 0.33 GHz, and appears 
to emit a series of bursts of typical duration ~ 10 min. 
The flux in a single burst is ~ 1 Jy and no emission is 
detected between the events to an upper limit of 75 mJy. 
The bursts seem to follow a regular pattern, with a time 
separation of ~ 77 min. The bursting source is unresolved 
and no counterparts are found thus far in other spectral 
bands. 

The quite high radio flux, combined with the limited 
spatial extent of the emitting region derived from the de- 
cay time of the bursts, implies a brightness temperature of 
Tb ~ 10 16 K placing the source at the GC distance. If this 
is the case, a coherent mechanism is required to power the 
observed radio emission (Hyman et al. 2005). 

A number of possible interpretations for the peculiar 
properties of GCRT J 1745-3009 have been examined by 
Hyman et al. (2005). Their main conclusion is that, al- 
though the option that this is a close-by radio-source (a 
very cold dwarf, or an extrasolar planet) can not be dis- 
missed, it appears much more likely that GCRT J 1745- 
3009 is indeed located close to the GC. Along this line, 
they consider both a radio-pulsar and a neutron star- 
neutron star binary as possible options. The case for a 
pulsar nature of GCRT J 1745-3009 is more thoroughly 
discussed in a recent paper by Zhu & Xu (2005). 

If one interprets the interburst time as due to orbital 
motion, according to the original suggestion by Hyman et 
al. (2005), the system may be similar to the double pul- 
sar J0737-3039 (Burgay et al. 2003; Lyne et al. 2004). In 
this Letter we further explore the possibility that GCRT 
J1745-3009 is a double neutron star system, following the 



proposal by Turolla & Treves (2004, hereafter TT04), that 
the continuous emission from J0737-3039 could be due to 
coherent radiation, produced at the shock formed by the 
interaction of the relativistic wind of the more energetic 
pulsar with the magnetosphere of the companion. 

2. THE MODEL 

We develop a simplified model for GCRT J1745-3009 
by assuming that: i) the source is at the distance of the 
GC, D — 8.5 kpc, ii) the interval between the bursts is 
an orbital period, P = 77 min, and hi) the binary sys- 
tem consists of two neutron stars (NS). Eight examples of 
NS+NS binaries in which one of the star is a radio pulsar 
are presently known. Among these J0737-3039 is unique 
since both its components are detected as pulsars. More- 
over J0737-3039 is the only NS+NS binary (P orb = 2.4 h, 
a ~ 9 x 10 10 cm, e ~ 0.09; Lyne et al. 2004), which exhibits 
a continuous radio emission (-Fi.4GHz ~ 5 mJy) about two 
times stronger than the emission from the two pulsars 
(Lyne et al. 2004). The fact that unpulsed radio emis- 
sion occurs in a double neutron star binary is of particular 
relevance in connection with GCRT J 1745-3009 and is at 
the basis of the scenario we are proposing. So far it has 
not been possible to establish if the continuous radio flux 
from J0737-3039 exhibits an orbital modulation. Were the 
continuous radio emission steady or weakly modulated, de- 
spite the alleged analogies, the two systems must differ in 
some respects, as discussed further on. 

The origin of the continuous emission from J0737-3039 
is still uncertain. The starting point is assumed to be 
the interaction of the wind of the most luminous pulsar A 
(P A = 0.02 s, E A ~ 6 x 10 33 ergs" 1 , B A ~ 6 x 10 9 G) 
with the magnetosphere of pulsar B (P B = 2.77 s, Eb ~ 
2xl0 30 ergs" 1 , B B ~ 10 12 G 4 ). A shock is produced when 
the energy density of the relativistic wind of A equals the 
magnetic energy density of B. Since Ea ^> Eb, the wind 
penetrates deep into the magnetosphere of B and the shock 



1 Department of Physics, University of Padova, via Marzolo 8, 35131 Padova, Italy; turolla@pd.infn.it 

2 Istituto Nazionale di Astrofisica, Osservatorio Astronomico di Cagliari, loc. Poggio dei Pini, Strada 54, 09012 Capoterra, Italy; pos- 
senti@ca.astro.it 

3 Department of Physics and Mathematics, Universita dell'Insubria, Via Valleggio 11, 22100, Como, Italy; treves@mib.infn.it 

4 Actually, the magnetic field of pulsar B could be somewhat smaller than this value, see Lyutikov (2004) 



1 



2 



Turolla, Possenti, & Treves 



is well within B's light cylinder radius at rg ~ 6 x 10 9 cm 
from B (Lyne et al. 2004; TT04). Under the reasonable 
assumption that the shock is the site of particle accelera- 
tion, TT04 suggested two possible scenarios for the con- 
tinuous emission. If the relativistic electrons originating in 
the shock fill a large volume (size s=s 10 13 -10 14 cm corre- 
sponding to w 1" at the estimated distance of J0737-3039, 
~ 0.6 kpc), optically thin synchrotron radiation is able to 
account for the observed flux. On the contrary in case 
the size of the emitting region is comparable to that of 
the shock itself, the only way to circumvent severe self- 
absorption is to consider coherent emission. 

In the picture we propose for GCRT J 1745-3009 the 
coherent emission is activated when the relativistic wind 
of one component (which we call A) penetrates the mag- 
netosphere of the other (B), in analogy with what likely 
occurs in J0737-3039. However, as already noticed and 
at variance with the case of J0737-3039, the strong radio 
emission has a finite duration, T ~ 10 min. A possible ex- 
planation is that the two NSs in GCRT J 1745-3009 arc in 
rather eccentric orbits, so that the shock is inside B's light 
cylinder only when the separation between the two stars 
drops below a limiting value. Clearly such a situation will 
occur only close to periastron, as illustrated in Figure 1. 

The orbit semi-major axis follows from Kepler's law 
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where Mtot — Ma + Mb ~ 2.6Mq is the total mass. 
The area swept by the radius in a time T is simply 
S = ira 2 \/l — e 2 (T/P), where e is the orbit eccentricity. 
From this it follows that the portion of the orbit centered 
at the periastron and travelled in a time T is limited by 
the values of the phase angle 4> and 2ir — <j>, where as a 
function of e is implicitly given by 
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The solutions of eq. (2) for different values of e and 
T/P = 10min/77min = 0.13 are shown in Figure 2 (full 
line). 

The distance of the shock from B, rs, is given by 
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where E>b = E>b{Rb){Rb / r s) 3 assuming a dipole field, 
Rb = 10 6 cm is the radius of star B, Bb(Rb) is the polar 
surface field, Ea is the wind luminosity emitted by A in 
a solid angle ^ttxa and r = a(l — e 2 )/(l — ecoscj)) is the 
orbital separation. 

The burst radio luminosity at v = 0.33 GHz is 

Lr ~ 4tt X rF^D 2 ~ 2.7 x 10 31 X i? ergs" 1 (4) 

where \R is the fraction of the total solid angle covered by 
the coherent radio beam. We take Lr to be proportional 
to the fraction of Ea intercepted by the shock 
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£ represents the efficiency for the production of (coher- 
ent) radio photons. Note that the radio luminosity which 
appears in the previous equation depends on the orbital 
phase through r and rs- To derive an estimate of Ea, we 
evaluate eqs. (3) and (5) at the periastron, where (f> = it, 
r = a(l - e)~3x 10 10 cm, and we take Lr(tt) = Lr, the 
latter as given by eq. (4). This results in 
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Inserting the previous expression for Ea into eq. (3) and 
imposing that the shock is at the light cylinder radius of 
B, rs = r c ,B = cPb/2tt ~ 4.8 x 10 9 Pb cm, results in a 
further relation between the phase <j) an d e 
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The request that rs = r c> s at a given point P guarantees 
that rs < r c ,s at all points along the orbit which are closer 
than P to periastron. 

In order for our model to work, conditions (2) and (7) 
must be satisfied together. Only in this case the shock is 
inside the light cylinder radius of B for the observed dura- 
tion of the burst T. The values of <j> solutions of eq. (7) are 
shown in Fig. 2 as a function of the eccentricity and for 
different values of Pb and Bb(Rb) = 10 12 G; all other pa- 
rameters have been taken equal to their reference values. 
As it can be seen from Fig. 2, the simultaneous fulfillment 
of both conditions is possible only in a quite restricted 
range of spin periods of star B, 0.31 s < Pb < 0.35 s. The 
corresponding values of the eccentricity are in the range 
~ 0.3-0.6. In an even tighter interval of periods, around 
Pb = 0.32 s, a second solution for <f) is present at higher 
values of the eccentricity. The allowed periods increase 
with B B ; however, for 5 x 10 11 G < B b (Rb) < 3 x 10 12 G 
they arc still in the 0.3-0.4 s range. Increasing (decreasing) 
the beaming factor \R (the efficiency £) results in shorter 
periods. The shock distance from star B, as derived from 
eq. (3), is shown in Figure 3 (left panel) for e = 0.4. By 
combining eqs. (3) and (5) it easy to see that LR(<j>) oc r§. 
The corresponding light curve is shown in Figure 3 (right 
panel); no radio emission occurs when star B phase angle 
is outside the range for which rs < r Ct B- In Figure 3 the 
orbital variation of rs and Lr is shown for the parameters 
typical of J0737-3039 for comparison. In this case it is al- 
ways rs < r c _B- The shock distance changes only slightly 
along the orbit because of the lower eccentricity and the 
orbital modulation of the unpulsed flux is ~ 10%. 

The wind luminosity of pulsar B can be estimated as 
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Although the ratio of the wind luminosities Ea/Eb is 
somewhat smaller than that inferred in J0737-3039, it is 
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still Eb <C Ea and the energy density of the wind of A ex- 
ceeds that of B by a factor ~ 10 at B light cylinder. This 
ensures consistency with our starting assumption about 
the shock forming as the result of the interaction of the 
wind of A with the magnetosphere of B. 

3. DISCUSSION 

The scenario outlined in the previous section shows that 
if one chooses for the model parameters values close to 
those introduced above, the main observational features 
of GCRT J1745-3009 can be successfully reproduced. Al- 
though the general picture we propose resembles that of 
J0737-3039, the two binary systems differ in more than 
one respect. In particular, the shock radiative efficiency 
needs to be higher in GCRT J1745-3009, - 1(T 3 as com- 
pared to - 1(T 4 in J0737-3039 (TT04). Lower values of 
£, in fact, result in smaller periods for pulsar B: this pro- 
duces a steep increase of Eb which soon would become 
~ Ea at B's light cylinder radius. Albeit coherence seems 
indeed to be required in both systems, no detailed model 
for coherent radio emission has been put forward so far. In 
TT04 some possibilities were suggested but the very model 
of coherent radio emission should be revisited possibly in 
the light of most recent studies on this topic (e.g. Fung & 
Kuijpers 2004; Zhang & Loeb 2004). 

flyman et al. (2005) reported that GCRT J1745-3009 
was not detected in two ~ 6 hr images obtained in 1996- 
1998 and in - 1 hr observations in 2002-2003. While no 
conclusive evidence can be drawn from the shorter obser- 
vations which have duration less than the orbital period, 
the lack of activity in 1996-1998 might pose a difficulty to 
the NS+NS binary model, as the same authors state. Here 
we only mention that an "on" -"off" state may be achieved 
if star A precesses. For \A ~ 0.1, the wind beam is <~ 70° 
wide and even a moderate precession angle can prevent 
the beam from A from intercepting the magnetosphere of 
B which subtends an angle < 10°. Interestingly, the char- 
acteristic time 1/fl ~ 116P 5 / 3 M- 2 / 3 (l - e 2 )s (D amour 
& Ruffini 1974) for changing the system geometry due to 
geodetic precession of the spin axis of pulsar A is ~ 3 yr, 
assuming equal masses (M — 1.3M Q ) for the two stars. 

The range of Pb and E>b{Rb) for which our model holds 
implies that neutron star B may be an active radio-pulsar. 
The required spin down luminosity of neutron star A (see 
eq. [5]) is also compatible with it being a canonical radio- 
pulsar with B A ~ 10 12 - 10 13 G and P A ~ 0.1 - 1 s or a 
mildly recycled pulsar with B A ~ 5 x 10 9 - 5 x 10 10 G and 
Pa ~ 10 — 20 ms. The latter possibility fits better with 
the standard scenario for the formation of a double neutron 



star binary (e.g. van den Heuvel & de Loore 1973), which 
predicts that the system should include a relatively young 
neutron star (the one resulting from the second supernova 
explosion) and a recycled pulsar (the first born neutron 
star), spun up to tens of millisecond period by the short 
phase of mass transfer from the evolving companion. 

Dctcctability of the pulsed signal from these two sources 
depends on the radio luminosity and favorable orientation 
of their beams. Given the position in galactic coordinates 
(I = 358.9°, b — —0.5°) and for a distance in the range 
6-12 kpc, the available models for the distribution of the 
ionized components in the interstellar medium (Taylor & 
Cordes 1993; Cordes & Lazio 2001) indicate that scatter- 
ing would prevent detection of any pulsating radio signal 
with P shorter than few seconds at the frequency of 330 
MHz. However, observations at frequencies higher than 1 
GHz would preserve the possibility of discovering a ~ 0.3 s 
pulsar, while observing at frequency > 2 GHz is required 
for detecting a ~ 10 ms pulsar. 

An X-ray luminosity coincident with the shock of Lx ~ 
10 32 ergs -1 should be present, assuming that ~ 10% of 
the wind luminosity intercepted by the magnetosphere is 
converted into X-rays. However, the expected flux is ex- 
tremely low, w 10 -4 ctss -1 for XMM-EPIC and Chandra- 
ACIS, assuming a power-law spectrum with index 2 and 
Nh ~ 10 22 cm~ 2 . Wc do not expect substantial emis- 
sion in other energy bands, as in the case of J0737-3039, 
although the younger NS may be associated with a super- 
nova remnant. 

Indirect support in favor of a double NS system for 
GCRT J1745-3009 could come from further radio obser- 
vations of J0737-3039. In fact, evidence that the continu- 
ous source is unresolved at the arcsec level would rule out 
the transparent scenario (TT04), leaving coherent emis- 
sion from the shock as the only viable option. If the pic- 
ture of GCRT J1745-3009 proposed here will be confirmed, 
one should consider that other similar systems may exist 
in the Galaxy. Their search among radio transients with 
programs analogous to that of Hyman et al. (2005), may 
be of the paramount importance for estimating the coales- 
cence rate on NS+NS systems, a vital issue in connection 
with the response of new generation gravitational wave 
detectors like VIRGO and LIGO. 
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Fig. 1. — A cartoon of the putative double NS system in GCRT J1745-3009. The orbit of star B relative to star A has e = 0.4. The dashed 
and full circles represent the position of B's light cylinder and of the shock, respectively (not to scale). Along the portion of the trajectory 
close to periastron and marked with a heavy line the shock is inside the light cylinder. The shaded areas show the wind of A (for illustrative 
purpose only). 




Fig. 2. — The solutions of eqs. (2) (full line) and (7) (dashed lines) as a function of the eccentricity. Each dashed curve is labelled by the 
corresponding value of the spin period Pb in seconds. 
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Fig. 3. — Left: The variation of the shock radius with the orbital phase. Actually, a shock can form only outside the hatched region, where 
rs < r c b- Right: Same for the unpulsed radio luminosity (in arbitrary units; see text for details). The dashed lines show the same quantities 
for J0737-3039. 



